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We investigate the effect of noise strength on the macroscopic ordering dynamics of systems with 
symmetric absorbing states. Using an explicit stochastic microscopic model, we present evidence for 
a phase transition in the coarsening dynamics, from an Ising-like to a voter-like behavior, as the noise 
strength is increased past a nontrivial critical value. By mapping to a thermal diffusion process, 
we argue that the transition arises due to locally-absorbing states being entered more readily in the 
high-noise regime, which in turn prevents surface tension from driving the ordering process. 
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Model systems with absorbing states have featured 
prominently in the quest for a fundamental understand- 
ing of far-from-equilibrium phase transitions and critical 
phenomena [1, 2 . An absorbing state is one that, once 
entered, cannot be exited. An excellent example is the 
extinction of a species. Further instances can be found 
in the dynamics of catalytic reactions [3 and of calcium 
channels in living cells [4 . 

A clear picture of how systems with multiple absorb- 
ing states behave is still emerging. These include cases 
where several competing species can drive each other to 
extinction. The static phase diagram for systems with 
two symmetric absorbing states [5]-[8] has recently been 
shown to have a much richer structure than those with a 
single absorbing state [1 , featuring in particular order- 
disorder transitions of several different types [6]. 

In this Letter, we turn to the dynamics of these sys- 
tems, focusing on how an absorbing state of global order 
is reached over time. We report the existence of a dynam- 
ical phase transition between two distinct macroscopic 
ordering regimes at a nontrivial value of the noise ampli- 
tude. More precisely, as the scale of local fluctuations in 
the dynamics is increased (whilst holding fixed all deter- 
ministic contributions) the surface tension that normally 
drives the coarsening process suddenly vanishes, leaving 
fluctuations as the only means to order. These two order- 
ing modes are characteristic of the Ising and voter models 
respectively [6 . The physical origin of this transition is 
elucidated by mapping the nonequilibrium stochastic dy- 
namics onto a thermal diffusion process, a technique pre- 
viously used for models with a single absorbing state [9] . 
We find that strong noise allows locally-absorbing states 
to be entered easily. Since the deterministic forces due to 
surface tension do not act in these states, further order 
can only be attained through interfacial fluctuations. 

This noise-induced dynamical transition is intriguing 
for several reasons. First, it has not been pre-empted 
in models with symmetric absorbing states: a renormal- 
ization group treatment [8 suggests that noise is an ir- 
relevant parameter as it is in the equilibrium theory of 
phase ordering kinetics [10]. Whilst specific microscopic 
models exhibit a wide range of phenomena [Q EJ |7l E], 



the relative strengths of the deterministic and stochastic 
contributions to the dynamics are often not independent. 
Here, we construct a model that avoids this shortcoming, 
and thus show that a change in noise strength alone can 
induce a dynamical phase transition. As we discuss in 
more detail below, our model is applicable to human lan- 
guage behaviour [TTl [12] in which noise strength relates 
to a memory decay rate [T2 . Thus, the spatial structure 
of linguistic diversity could in principle be nontrivially 
affected by how quickly humans forget. Finally, our find- 
ings bear on basic questions about nonequilibrium uni- 
versality classes, a point we return to in the conclusion. 

Our starting point is the Langevin equation proposed 
in Ref. [6 for the general class of systems with two ab- 
sorbing states. It reads 

dt(^ = {a^-b^^){l - (/)2)+W'(/) + aVl-(/)2 rj{t) , (1) 

where (j) G [—1,1] is a continuous coarse-grained field, 
usually the magnetization. The first term is a determin- 
istic force that derives from a potential; the second leads 
to smoothing of interfaces through surface tension (see 
e.g. [10 ); and 77 is a Gaussian white noise with zero mean 
and unit variance. A key feature of this equation is that 
all terms vanish in (open) regions where = dzl. Such 
regions are locally absorbing: changes in may then only 
take place at domain boundaries. 

The static phase behavior is found to depend on the 
shape of the potential, i.e., the values of a and 6 in ([T]) 
[6 . Here we are interested in what happens when the 
potential is held fixed, but the noise strength a is varied. 
As an alternative to direct numerical integration of the 
Langevin equation ([T]), the method used in [6 , we con- 
struct a stochastic microscopic model that is described 
by this Langevin equation at the mesoscopic scale. 

We first define the dynamics of the microscopic model, 
and then show that it is described by a spatially-discrete 
version of ([T]). It is defined on an L x L square lattice, 
each site of which hosts N binary spins. In each update, a 
spin is chosen at random and replaced with a copy from 
its neighborhood. To realize the diffusion term in ([T]), 
there is a probability q = ^ that the copy is taken from 
one of the z = 4 neighboring sites. Otherwise, the copy is 
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FIG. 1: Dynamics of the microscopic model. Each lattice 
site (shown as a region bounded by a dashed line) contains 
N spins. In each update, a randomly-chosen spin is replaced 
with a copy taken from one of the z = 4 nearest-neighbor 
sites, each with probability q/z; an up-spin from the same 
site with probability (1 — q)p or a down-spin from the same 
site with probability (1 — ^)(1 — p)- 



taken from the same site. The potential term in ([T]) can 
be interpreted as a systematic bias towards copying one 
of the two spin states when copying from the same site. 
Specifically, if 6 = 0, one should choose an up-spin with 
probability p = ^(l + ^i) + ^0i(l — 0f), where is the 
local magnetization. That is, = {n^^i — ni^i)/N ^ where 
n^^i and n^^i are the numbers of up and down spins on 
site i. These update rules are illustrated in Fig. [T] 

One concrete application of this model is in linguis- 
tics p^. Each site represents a speaker, and each spin 
represents a previously-heard utterance stored in mem- 
ory. The two spin states relate to two different ways of 
saying the same thing (e.g., phonetic realizations of a 
vowel). The local bias, in which the majority spin state 
is favored, then models an experimentally-observed ten- 
dency for language users to overproduce the most fre- 
quent variant [11]. The key point is that the number of 
spins per site, sets both the noise strength (see below) 
and the lifetime of an utterance in memory. An appro- 
priate choice for this lifetime was a central consideration 
in an analysis of the New Zealand English dialect [12 . 

The stochastic equation of motion for the system is 
obtained by calculating moments of the change S<pi in 
the local magnetization after one update [EllM]. All of 
these jump moments vanish in the limit ^ oo as some 
power of 1/N. The limits limAr^oo A^^(^^i) = A({0}) 
and limAT^oo ^^(^^f ) = ^iiW}) finite, in which 

Ai = a(/.i(l - 4>f) + hJ2i4>j - 'Pi) , Bi = 2(1 - 4>f) , 

j 

(2) 

and the sum is over nearest neighbors of site i. All other 
jump moments vanish at least as fast as A Fokker- 

Planck equation for the model is then obtained from a 
Kramers-Moyal expansion [13 . We define a time step 
as 5t = and take the limit A' ^ oo to obtain a 

continuous time limit. The Fokker-Planck equation 
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ning stochastic simulations of the microscopic dynamics 
we set A^ = 100 to allow for practicable run times. 

This Fokker-Planck equation can be written equiva- 
lent ly (and thus also exactly) as a Langevin equation 
under the Ito prescription fT3|. We have 



dt^, = h(rcl,,{l-^'i)^yU,-cP, 
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exactly describes the stochastic dynamics of our micro- 
scopic model in this continuous-time limit. When run- 



which is a spatially- discrete version of Eq. ([T]) with a = 
r/i, b = 0^ D = h5'^^ where S is the lattice spacing, and 
<j = 1. Written in this way it is clear that if h is varied 
while holding r constant then the strength of the noise, 
relative to the deterministic term, can be varied without 
changing the shape of the potential. We thus study the 
effect of noise strength in our Monte Carlo simulations 
of the microscopic dynamics by varying h. In all these 
simulations we take r = 3, since we have found this to 
be large enough to provide the metastable interfaces be- 
tween regions of positive and negative (j) that are required 
if domain coarsening is to take place [18]. 

We now examine the dynamics in the high (small h) 
and low (large h) noise regimes. On a finite system, and 
with a > 0, 6 = in ([T]), the system eventually ends up 
in one of the globally absorbing states (i.e., all <pi = 1 or 
all = —1) In the large- regime, we anticipate a 
domain growth driven by surface tension. This is because 
the deterministic limit h ^ oo corresponds to the time- 
dependent Ginzburg-Landau equation, obtained from the 
Landau free energy functional for the Ising model with 
a non-conserved order parameter When the noise 

amplitude is small, we do not expect the distinction be- 
tween additive and multiplicative noise to be important, 
and thus the dynamics should coincide with model A for 
the Ising model. On the other hand, it can be shown that 
in the limit ^ 0, the purely fluctuation-driven dynam- 
ics of the voter model are formally recovered [15]. This 
suggests the possibility of at least a crossover (if not a 
transition) from Ising to voter coarsening as h is reduced. 

One way to identify the coarsening regime is to exam- 
ine the density of interfaces as a function of time [5] [6] . 
This is defined as p{t) = 4j;2 j) ~ ^i^j) ^tie 
square lattice in two dimensions, in which the sum is 
over distinct nearest-neighbor pairs. In the Ising model, 
p{t) - [To], while in the voter model, p{t) - l/ln(t) 

in two dimensions fT6^. We measure p{t) by averaging 
over multiple realizations of the microscopic stochastic 
dynamics described above. The results, shown in Fig. [2| 
are suggestive of algebraic coarsening taking place at 
large /i, and logarithmic coarsening at low h. We remark 
that the deviation from the expected law seen at 

late times for large is a finite-size effect caused by a 
domain coalescing with a periodic image of itself. 

We now examine whether this shift from Ising to voter 
coarsening is a sharp transition or a smooth crossover, 
and if the former, what the mechanism for a transi- 
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FIG. 2: (Color online) Interfacial density p as a function of 
time t in Monte Carlo simulations of the microscopic model. 
Left: ^/p{t) is plotted on logarithmic axes, such that a 
straight line indicates algebraic coarsening. For large h, and 
at times before the onset of finite-size effects, the gradient 
is consistent with the Ising model value of | (dashed line). 
Right: [p{t) lnt]~^ plotted on linear axes. A constant asymp- 
tote, seen for small /i, indicates logarithmic coarsening. 



tion would be. The multiplicative noise in the Langevin 
equation Q makes interpretation of the stochastic dy- 
namics difficult as it stands. Our physical intuition is 
strongest when the noise is additive, and the dynamics 
can be viewed in terms of diffusion in a potential. To 
this end, we transform the local magnetization 0^ to a 
variable Oi such that the Fokker-Planck equation analo- 
gous to (|3| has a diffusion term Bi{{0}) that is indepen- 
dent of the coordinates 0^. The appropriate transforma- 
tion is Oi = sin^^ (see p!7| for details). The resulting 
Fokker-Planck equation corresponds uniquely to the set 
of Langevin equations 0^ = —V^{6i) + r^, one for each 
site i, and where r] is the usual thermal (Gaussian white) 
noise. The potential Vd has the form 
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where m is the mean magnetization of the z sites that 
are neighbors of site i. Note the separate Langevin 
equations are coupled because m depends on the values 
of at neighboring sites. 

The subscript D here is used to emphasize the cru- 
cial difference between this potential, felt by a thermal 
diffusion process, and that whose derivative gives the de- 
terministic term in Q . The key points are that the shape 
of Vd gives insight into the stability of a microscopic con- 
figuration under the stochastic dynamics, and that this 
shape changes with the noise strength in a nontrivial 
way. For our purposes it suffices to look at the behavior 
near the boundaries, from which we deduce three distinct 
qualitative shapes that this diffusion potential may take. 

As approaches either boundary point ^ = ±| at 
some fixed local magnetization m, Vd diverges logarith- 
mically. Depending on the values of h and m, the diver- 
gence can be either towards +oo or — oo. Considering the 
case m < 0, we find that the divergence is towards +oo 
at the right boundary when h > h- = [22^(1 — m)]~^. 
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FIG. 3: (Color online) Shape of the diffusion potential at 
increasing h. Here, h± = [22;(1 zbm)]~^, z = 4 and fh = —0.5. 



and towards — oo for smaller values of h. Similarly, at 
the left boundary, the divergence is towards +oo when 
h > hj^ = [2z{l -\- fh)]~^ . The three possible combi- 
nations of boundary divergence are as shown in Fig. [3] 
Analogous shapes are found for the case m > by using 
the symmetry Oi —Oi and fh ^ —fh. 

We contend that this changing shape of the diffusion 
potential is of significance for the coarsening dynamics 
for the following reason. For any microscopic spin con- 
figuration, we always have that h- > j^. Therefore, for 

< there are always minima of the diffusion poten- 
tial at the boundaries, and a single maximum in between 
(since there are at most two extrema in the interior re- 
gion). Thus the diffusion of is at any instant biased 
towards a locally absorbing state. Once a locally absorb- 
ing state is reached, the deterministic force in Eq. Q 
vanishes. Then, all that remains is the Langevin equa- 
tion for the voter model (the form of which is given e.g. in 
) . This leads us to suggest that voter coarsening should 
be seen over some finite range of /i, at least up to h = 
Since the deterministic limit of Q and model A for tfie 
Ising model are equivalent, there is the possibility of a 
transition to Ising-type coarsening at some nonzero value 
of /i, perhaps in the vicinity oi h = ^. 

This picture is confirmed if we can identify a point 
at which the ordering dynamics changes from the voter- 
type to the Ising-type from Monte Carlo simulation data. 
Here we have found a droplet experiment^ similar to that 
described in [5j, most useful. This entails an initial con- 
dition in which all spins within a radius Rq of the origin 
are up, and the rest are down. When coarsening under 
surface tension, the Allen-Cahn theory of phase ordering 
predicts that the droplet will shrink, its area decreasing 
at a constant rate proportional to the diffusion constant 
D Since D (x h in our model, we expect the mag- 

netization density m{t) = j2 decrease at a rate 

c(x.h/L'^ through surface tension. On the other hand, in 
the voter model, the rate of change of the magnetization, 
averaged over multiple stochastic realizations, is zero. 

To each stochastic realization of m(t) we fit a linear 
function. An estimate of c, and an error, can be com- 
puted from the mean and standard deviation of the mea- 
sured gradients. We can compare results from different 
system sizes L by taking the initial droplet radius Rq 
proportional to L, and by plotting cLF' against h. These 
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FIG. 4: (Color online) Droplet shrinking rate, c, as a function 
of inverse noise strength /i at a range of system sizes L. Each 
data point is obtained from a sample of 100 simulation runs. 



results are shown in Fig. [4j The expected linear increase 
of c with h is observed. What is interesting is that the 
intercept is not at = 0, but at some positive value of /i, 
below which c is consistent with zero. This suggests that 
there is indeed a transition between voter and Ising coars- 
ening dynamics at some h = > 0. A least-squares 
linear fit to the data yields an estimate of /i* ^ 0.059, 
which is close to the value ^ = 0.0625 suggested by the 
analysis of the diffusion potential. It is possible that a 
nonzero transition point could be a finite-size effect. The 
fact that the data show no systematic shift of /i* towards 
zero as L is increased appears to rule this out. 

To summarize, we have presented evidence for a noise- 
induced transition in the ordering dynamics in systems 
with two symmetric absorbing states and discrete space. 
The origin of this transition is that, once the noise is 
sufficiently strong, /oca//^-absorbing states are entered 
with ease. This eliminates the term in Q that leads 
to surface-tension driven coarsening. Then, fluctuation- 
driven voter-type coarsening is the only possible way the 
system can order. We believe discrete space is essential 
for the operation of this mechanism, so that a locally- 
absorbing state = ±1) can be reached more quickly 
than the diffusion process connecting neighboring sites. 
By contrast, in continuous space, any finite transition 
point /i* > collapses onto a zero diffusion constant, 
D = 0. Therefore, one would generically expect Ising- 
type ordering on the continuum. This could explain why 
a previous renormalization group treatment of ([T]) did 
not reveal the noise strength as a relevant quantity [8] . 

However, the fact the combination of discrete space 
and non-thermal noise leads to a further phase transition 
is in itself interesting. It would be worthwhile to try and 
understand this transition more rigorously than through 
the heuristic diffusion potential picture we have described 
here. The fact that noise strength is related to memory 
lifetime in the application of our microscopic model to 
language change also warrants further investigation, in 



particular, whether the transition is also seen on social 
network structures and whether the resulting patterns of 
diversity can be distinguished empirically. 

We conclude by returning to a question investigated 
by 0: that is, which qualitative features, such as sym- 
metries and conservation laws, determine the dynami- 
cal universality class of systems with multiple absorbing 
states. We have seen that a set of models with the same 
qualitative properties, but different noise strengths, can 
have different ordering dynamics. In particular, strong 
noise leads to the magnetization conservation character- 
istic of the voter model being an emergent consequence 
of the stochastic dynamics. To develop a better general 
understanding of nonequilibrium phase transitions and 
critical phenomena, one may need to determine whether 
other cases exist in which such conservation laws emerge. 
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